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. Abstract. Equation of motion of an uncharged arbitrarily shaped dust particle 

under the effects of (stellar) electromagnetic radiation and thermal emission is 
derived. The resulting relativistically covariant equation of motion is expressed in 
terms of standard optical parameters. Relations between energy and mass of the 
incoming and outgoing radiation are obtained, together with relations between 



radiation energy and mass of the particle. The role of the diffraction nicely fits 
• ■ the relativistic formulation of the momentum of the outgoing radiation. 

Covariant formulations yield several simple consequences. It is shown that 



the frequently used statement "energy equals 7 x mass x c 2 " or "energy equals 
mass x c 2 " are not correct, in general; "rest energy equals mass x c 2 " holds. 
The inequality < Q pr ,i/Q'ext < 2 is a simple relativistic consequence for the 



Poynting- Robertson (P-R) effect (Q' ext and Q' 1 are dimensionless efficiency fac- 
tors for the extinction and radial direction of the radiation pressure, integrated 
over stellar spectrum). The condition for the P-R effect is p' D = ( 1 — Q' pr i/Q'ewt 
) p'i, where p\ and p' a axe incoming and outgoing radiation momenta (per unit 
time) measured in the proper frame of reference of the particle. The case of "per- 
fectly absorbing spherical dust particle" , within geometrical optics approximation, 
corresponds to the condition p' = 0.5 p'i. The obtained results enable understand- 
ing of the physics of the P-R effect and evaluate the published explanations. As 
for arbitrarily shaped dust particle, the general condition < C' pr j / C' ext < 2 
/ (1 + 5Zj=2 ^pr, j /C'pr,!) holds for cross sections of extinction and radiation 
pressure components, if thermal emission force is of negligible importance. The 
condition can add a new information to the results obtained from observations, 
measurements and numerical calculations of the optical properties of the particle. 



Key words, cosmic dust, electromagnetic radiation, thermal emission, relativity 
theory, equation of motion 
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1. Introduction 



The Poynting- Robertson effect is used in modelling orbital evolution of dust grains un- 
der the action of electromagnetic radiation (of the central star), for many decades (e. g., 
Poynting 1903, Robertson 1937, Wyatt and Whipple 1950, Dohnanyi 1978, Kapisinsky 
1984, Jackson and Zook 1989, Leinert and Grim 1990, Gustafson 1994, Dermott et al. 
1994, Reach et al. 1995, Murray and Dermott 1999, Woolfson 2000, Danby 2003, Quinn 
2005, Harwit 2006, Gran 2007, Sykes 2007, Kriigel 2008). It is assumed that particle 
with spherically distributed material can be used as a correct approximation to reality. 
However, real nonspherical particles interact with electromagnetic radiation in a compli- 
cated manner and particles of various optical properties exist (e. g., Mishchenko et al. 
2002). Thus, it is essential to have an equation of motion sufficiently general to cover a 
wide range of optical parameters, not just the limited cases standardly investigated. The 
first presentations of the equation of motion for arbitrarily shaped particle were given by 
Klacka (1994), Klacka and Kocifaj (1994), Klacka (2000), Kocifaj et al. (2000), Klacka 
and Kocifaj (2001). Relativistically covariant derivation can be found in Klacka (2000) 
and Klacka (2004), where also thermal emission of the particle is taken into account - 
Mishchenko (2001) has formulated the radiation pressure on arbitrarily shaped particles 
arising from an anisotropy of thermal emission. Krauss and Wurm (2004) have published 
the first experimental evidence that nonspherical dust grains behave in a different way 
than the spherical ones. 

We derive equation of motion for arbitrarily shaped particle in a physically correct 
way and relativistic covariant formulations are presented. A new form for the validity of 
the Poynting-Robertson (P-R) effect is found: the P-R effect holds only when the total 
momentum of the outgoing radiation per unit time is colinear with the incident radiation, 
in the proper/rest frame of reference of the particle, and, a new physical condition for 
the P-R effect is found, as for the dimensionless efficiency factors for radiation pressure 
and extinction (optical properties) of the particle. Moreover, relations between energy 
and mass of the incoming and outgoing radiation are presented, both for arbitrarily 
shaped particle and also, as a special case, for spherical particle for which the P-R effect 
holds. The relations between energy and mass for the incoming and outgoing radiation, 
together with relations between radiative energy and mass of the particle, present a 
practical contribution to the discussion on relativistic energy-mass relation (see, e. g., 
Ugarov 1969, Sees. 57-58, pp. 246-252; Okun 1989a, 1989b; Taylor and Wheeler 1992, 
Sec. 8.4, pp. 228-233; Einstein 1999). The relations are consistent with the consideration 
of the diffracted light within geometrical optics approximation. The results of the paper, 
including detail discussion also with respect to the statements presented in the original 
papers by Poynting (1903) and Robertson (1937), enable to understand physics of the P-R 
effect, which is frequently used in studies of orbital evolution of cosmic dust particles. As 
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for the arbitrarily shaped dust particles, a condition yielding some information on cross 
sections, unknown from experiments, observations or theoretical/numerical solutions is 
obtained as a simple relativistic consequence of the covariant form of the four-momentum 
of the outgoing radiation. 

We begin by reviewing in Sec. 2 and 3 the basic physical processes in proper and 
stationary frames, for an arbitrarily shaped particle. The physical equations use cross 
sections for all relevant phenomena (extinction, scattering, absorption, radiation pres- 
sure components) and the equation of motion of the particle for the case of the irradiated 
particle is discussed at the end of Sec. 3 (Eqs. 40-41). The following Sec. 4 deals with the 
covariant formulation for the incoming and outgoing radiation, which are used in Sec. 5 
for finding energy-mass relation for the radiation. Sec. 6 discusses the relation between 
the change of the particle's mass, separately for the incoming and outgoing radiation. 
Application of the results from the previous sections are discussed, for special particle 
shapes, in Sec. 7: (i) it is shown that a simple relations holds for dimensionless efficiency 
factors for radiation pressure and extinction for the case of spherical particle (Eq. 79), 
including presentation of the value for perfectly absorbing (or perfectly reflecting) sphere 
(Eq. 104) which differs from the conventional statement coming back to Poynting and 
Robertson; (ii) an understanding of the action of the radiation on motion of spherical 
particles is discussed in detail in Sec. 7.4; (iii) an analogy between mechanics and elec- 
tromagnetism for the case of planar surface is discussed in Sec. 7.5. Results of Sees. 2-7 
are used in Sec. 8, which presents considerations on the equivalence principle of mass 
and energy for absorbing planar surfaces. A special Sec. 9 derives a simple condition for 
cross sections for radiation pressure components and extinction, using results of Sec. 5. 
The most important results are shortly summarized in Sec. 10. 

2. Proper reference frame of the particle - stationary particle 

The term "stationary particle" will denote a particle which does not move in a given 
inertial frame of reference. Primed quantities will denote quantities measured in the 
proper reference frame of the particle - rest frame of the particle. 

The flux density of photons scattered into an elementary solid angle dfl' = 
sin#' d6' d(f>' is proportional to p'(9',(/)') dQ! , where p'(6',(f>') is the "phase function". 
The phase function depends on orientation of the particle with respect to the direction 
of the incident radiation and on the particle characteristics; angles 6', <j>' correspond to 
the direction (and orientation) of travel of the scattered radiation, 9' is the polar angle 
which vanishes for propagation along the unit vector e[ of the incident radiation. The 
phase function fulfills the normalisation condition 




(1) 
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The momentum of the incident beam of photons which is lost in the process of inter- 
action with the particle is proportional to the cross section C' ext (extinction). The part 
proportional to C' abs (absorption) is emitted in the form of thermal radiation and the 
part proportional to C' ext — C' abs = C' sca is scattered. The differential scattering cross 
section dC' sca /dQ,' = C' sca p'(0', <j>') depends on the polarization state of the incident light 
as well as on the incidence and scattering directions (e. g., Mishchenko et al. 2002). 

The momentum (per unit time) of the scattered photons into an elementary solid 
angle dCl' is 

dp' sea = \S' C' sca p'(6', <jj) K' dn' , (2) 
where the unit vector in the direction of scattering is 

K' = cos9' e[ + sine' cose/)' e' 2 + sinfl' sin 4/ e' 3 . (3) 

S' is the flux density of radiation energy (energy flow through unit area perpendicular to 
the ray per unit time). The system of unit vectors used on the RHS of the last equation 
forms an orthogonal basis. The total momentum (per unit time) of the scattered photons 
is 

p' sca = - S' C' sca f p'{6',<t)') K' dSH . (4) 

c Jin 

The momentum (per unit time) obtained by the particle due to the interaction with 
radiation - radiation force acting on the particle - is 

% = \S' \c' ext e[ C' sca £ p'{6',4>') K> dO'J + F' e (T>) , (5) 

where the emission component of the radiation force acting on the particle of absolute 
temperature T" is (Mishchenko et al. 2002, pp. 63-66) 



K(T') = ~- r dJ [ r' K' e (r',T',J) dr' 

C Jo J4tt 



(6) 



The unit vector r = r' jr' is given by position vector r' of the observation point with 
origin inside the particle (the emitted radiation in the far-field zone of the particle prop- 
agates in the radial direction, i. c., along the unit vector f'), u>' is (angular) frequency of 
radiation, 



K' e {r',T',J)=I' b {T',u>) 



K' u (r',u')- [ Z' n (r',r",u;')dr" 
Jin 



(7) 



I in 

where K' u is the (1,1) element of the particle extinction matrix, Z' n is the (1,1) element 
of the phase matrix and the Planck blackbody energy distribution is given by the well- 
known relation 

Thermal emission has to be included in the total interaction of the particle with elec- 
tromagnetic radiation: if the particle's absolute temperature is above zero, it can emit 
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as well as scatter and absorb electromagnetic radiation. The particle is assumed to be 
isothermal, in Eqs. (6)- (8). 

Equation (5) can be rewritten to the form 

^ = \ S ' {[CLt ~ <cos6'> C" sca ] ei + 
[ - < sin 6' cos 4> > C' sca ] e' 2 + 

3 

[ - < sin^' sin0' > C' sca ] e' 3 } + ]T F' e - e> , (9) 

3=1 

where < x' >= J i7i x' p'{8', <fi') dVt' and F' e ■ = F' e (T') ■ e'y As for the energy, we assume 
that it is conserved: the energy (per unit time) of the incoming radiation E[, equals to 
the energy (per unit time) of the outgoing radiation (after interaction with the particle) 
E' Q . We will use the fact that time t' = r, where r is proper time. 

Summarizing important equations, we can write them in a short form 

3=1 

where C' prA = C' ext - < cos 9' > C' sca , C' pr 2 = - < sin#' cos0' > C' sca , C' pl%3 = 
— < sm8' sin<p' > C' sca are cross sections for radiation pressure. We have added an 
assumption of equilibrium state when the particle's mass does not change. 



2.1. Summary of the important equations 

Using the text concerning energy below Eq. (9) and the last Eq. (10), we may describe the 
total process of interaction in the form of the following equations (energies and momenta 
per unit time): 

E' = E[ = S' C" ext , 

^ext / V ^ext ^ext / c J=1 



V\ = f ei . (11) 

The index "i" represents the incoming (incident) radiation, beam of photons, the index 
"o" represents the outgoing radiation. 

The changes of energy and momentum of the particle due to the interaction with 
electromagnetic radiation are 

P'i- P'o- (12) 



dr 
dp' 



dr 

3. Stationary frame of reference 

By the term "stationary frame of reference" (laboratory frame) we mean a frame of 
reference in which particle moves with a velocity vector v = v(t). The physical quantities 
measured in the stationary frame of reference will be denoted by unprimed symbols. 
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Our aim is to derive equation of motion for the particle in the stationary frame of 
reference. We will use the fact that we know this equation in the proper frame of reference 
- see Eqs. (11) and (12). 

If we have a four- vector = (A° , A) , where A is its time component and A is its 
spatial component, generalized special Lorentz transformation yields 

A' = 7 (A - v-A/c) , 

A' = A + [( 7 - 1) v ■ A/v 2 - 7 A°/c] v , (13) 
with inverse 

A" = 7 (A + v ■ A'/c) , 

A = A' + [(7 - 1) v ■ A'/v 2 + 7 A'°/c] v , (14) 



where 7 = 1 / \J\ — v 2 /c 2 . 

As for four-vectors we immediately introduce four-momentum: 

p» = (p ,p)^(E/c,p) . (15) 
3.1. Incoming radiation 

Applying Eqs. (14) and (15) to quantity {E^/c^p'A (four-momentum per unit time 
proper time is a scalar quantity) and taking into account also Eqs. (11), we can write 

Ei = E'i 7 (1 + v ■ e'Jc) , 
E' 

Pi = -± {e[ + [(7 - 1) v ■ e[/v 2 + 7 /c] v} . (16) 
Using the fact that = (h vj 'c,h v e\/c) for photons, we have 

v 1 = V W\ , 

e'i = — {ei + [(7 - 1) v- ei /v 2 - 7/c] v} , (17) 

Wl 

where 

Wi = 7 (1 — v ■ ei/c) . (18) 

Inserting the second of Eqs. (17) into Eq. (16), one obtains 
Ei = (1/ Wl ) El , 

Pi - (1M) (E'Jc) ei . (19) 

We have four-vector p^ — (Ei/c,Pi) = (l,e{)Ei/c = (1/wi, ei/wi) x Wi Ei/c = 
Wl Ei/c. 

For monochromatic radiation the flux density of radiation energy becomes 
S' = n' hv' c ; S = nhu c, (20) 
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where n and n' are concentrations of photons (photon number densities) in the corre- 
sponding reference frames. We also have continuity equation 

^J" = °. f = (cn,cn ei ), (21) 
with current density j^. Application of Eq. (13) then yields 

ri = w 1 n . (22) 
Using Eqs. (17), (20) and (22) we finally obtain 

S' = wj S . (23) 
Eqs. (11), (19) and (23) then together give E t — w\ S C' ext , pi — w\ S C' ext ei/c. 

3.2. Outgoing radiation 

The situation is analogous to that of the preceding subsection. It is only a little more 
algebraically complicated, since radiation may also spread out in directions given by unit 
vectors e2, e 3 . We need transformations e'j — > ej, j = 2, 3. The vectors e 2 and e' 3 can be 
used to describe directions of propagation of radiation scattered by the particle. Thus, 
aberration of light also exists for each of these unit vectors. The relations between e' 2 
and e 2 , e' 3 and e 3 , are analogous to that presented by the second of Eq. (17): 



e 



1 -— {ej + [( 7 - 1) v-ej/v 2 - j/c] v} , j = 1,2,3, (24) 



W 



3 

where 

u>j = 7(l - vej/c) , .7 = 1,2,3. (25) 

It is worth mentioning that vectors {e'pj — 1,2,3} form an orthonormal set of vectors, 
and, unit vectors {ej;j = 1,2,3} are not orthogonal unit vectors. 

Applying Eqs. (14) and (15) to the quantity (E' /c,p' ) (four-momentum per unit 
time - proper time is a scalar quantity), we can write 



E = 7 (E' Q + v ■ p' ) 



Po = p' + 



v . (26) 



Using also p\ = E[ e'Jc and Eqs. (11), (24), (26), 



E = -^w lEtl + 1 - Ei 

^ext \ u ext / 

V C ext C ext ) V C ext w 2 O ext W 3 

3 

' C 



p° = ( 1 - ~rr- J — e i + 77, ^- 7 « 
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ECp wi Ei ( ej \ 



J=2 
3 



3 = 1 



3.3. Equation of motion 

In analogy with Eqs. (12), we have for the changes of energy and momentum of the 

particle due to the interaction with electromagnetic radiation 

dE 

—, — — Ei — E , 
d t 

ir = Pi - Po- (28) 

a r 

Putting Eqs. (27) into Eqs. (28), using also pt = ( Ei / c) e\, one easily obtains 



i=i 



Eq. (29) may be rewritten in terms of four- vectors: 
where p M is four-vector of the particle of mass m 

pf* = mu f " , (31) 
four- vector of the world- velocity of the particle is 

= (7 c,7 v) . (32) 
We have also other four-vectors 

6^ = (1/ Wj , ej / Wj ) , . 7 = 1,2,3. (33) 

It can be easily verified that: 

i) the quantity w Ei is a scalar quantity - see first of Eqs. (19); 

ii) Eq. (30) reduces to Eq. (10) for the case of proper inertial frame of reference of the 
particle; 

iii) Eq. (30) yields d m/d r = 0. 
We introduce 



= l/wj wl + v ■ e'j I c , 
= e j I Wj ~ e'j + v I c , 
3 = 1, 2, 3 



(34) 
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for the purpose of practical calculations. Physics of these relations corresponds to aber- 
ration of light. In general case, b'^ = (l,e^), b^ = A l/ ^b' j 1 ', j e {1,2,3}, where K v » 
represents general Lorentz transformation. Eq. (34) is correct to the first order in v/c, 
as it is presented. 

We have derived an equation of motion for real dust particle under the action of 
electromagnetic radiation (including thermal emission) . It is supposed that the equation 
of motion is represented by Eqs. (11) and (12) in the proper frame of reference of the 
particle. The final covariant form is represented by Eq. (30), or using Ei — w\ S C' ext 
(see Eqs. (11), (19) and (23)), 

j=i 

To first order in v/c, Eqs. (34)- (35) yield 
d i) S ^ 

= — E C ' P rJ K 1 - 2 v ■ e i/ c + v ■ e o/c) ej v/c] + 



i t K, K 1 + v ~ir) "i ' 



m 

3 = 1 



C / C- 

ej = (1 - v-e'j/c) e'j + v/c , j = 1,2,3. (36) 

It is worth mentioning to stress that the values of radiation pressure cross sections C' pr 3 , 
j = 1, 2, 3, depend on particle's orientation with respect to the incident radiation - their 
values are time dependent. 

It can be verified that Eq. (35) (or Eq. 36 within the accuracy to the first order in v/c) 
yields as special cases the situations discussed in Einstein (1905) and Robertson (1937). 
However, Eqs. (11) and (27) are not consistent with the statement presented by Poynting 
(1903), Robertson (1937), Wyatt and Whipple (1950) and others about the process of 
reemission of perfectly absorbing spherical dust particle (explanation is presented in Sees. 
7.3.3 and 7.3.4). 

3.4. Continuous distribution of density flux of energy 

For a continuous frequency distribution of density flux of energy, we can write 
dp' 



d T 

3 

S E C 'vr, + Ki) ■ ( 37 ) 

Taking into account that concentration of photons fulfills n' = w\ n (Eq. (22)) and 

that Doppler effect yields v' = w\ v (Eq. (17)), we have dn'/dv' = dn/dv. Lorentz 
transformation finally yields 



dp^ 
~d 



^ f w\ f 00 , dn . . , 1.1 
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= E(^ ^ + ^ ^) (c 6? - . (38) 

As a consequence, dm/dr — (this corresponds to the condition dE'/dr = 0). As for 
the accuracy to the first order in v/c, equation of motion of the type Eq. (36) holds if 
the substitution C j — > C' j, j = 1, 2, 3 is done: 

d i) S 

77 = — E ^ K 1 - 2 ^ • e i/ c + « • *i/c) ^ - t;/c] + 

i e [(' + ^) « - g ■ 

ej = (1 - v-e'/c)e' + v/c , j = 1,2,3. (39) 



The values of radiation pressure cross sections C' pr j, j = 1 to 3, depend on particle's 
orientation with respect to the incident radiation. 

If the particle is not irradiated for a time interval, then outgoing radiation due to the 
thermal emission is given by four- vector ( E' a /c, p' a ), where E' a is energy per unit time 
and p' Q = — F' e = — X^=i F'e j e 'j i s momentum per unit time. Using transformations 
represented by Eq. (14), or Eq. (26), and using also Eq. (24) together with the fact that 
the four- force acting on the particle is ( d p 11 / d r ) e = — p£, one obtains 

instead of Eq. (38) and E' Q is energy per unit time, which is lost due to the thermal 
emission. Eq. (40) yields 



= 1 V F' (V - ^ 

\&t L m f-f e ' j \ j c 



3= 

i • (41) 



- 



drj e 

Mass of the particle decreases due to the thermal emission, alone. 



4. Incoming and outgoing radiation - covariant formulation 

We have treated the process of interaction between the electromagnetic radiation and 
a dust grain in two ways: we have considered the incoming/incident radiation and the 
outgoing radiation. As we are interested in the relations between the corresponding part 
of radiation and it's mass, we formulate covariant forms for the two parts of radiation, 
in this section. 
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4.1. Momentum-four vector for incoming radiation 

On the basis of Eq. (11) and Sec. 3.1, mainly Eqs. (18)-(23), and also using the results 
of Sec. 3.4, we can write four-momentum of the incoming (incident) radiation 

dpi = dr , (42) 

where also Eqs. (25) (or 18) and (33) have to be used. 

4.2. Momentum-four vector for outgoing radiation 

On the basis of Eqs. (23), (25), (27), (32), (33) and the results of Sec. 3.4, we can write 
for the energy and momentum of the outgoing radiation 



X = ^ S ° Lt . (43) 

c 

Eq. (43) can be written in a short relativistically covariant form 



X = 



c 

- (^f , Pout J (44) 
and Eqs. (25), (32) and (33) can be used. 

4.3. Relation between outgoing and incoming radiation 

Eq. (44) reads 



Eqs. (42) and (45) imcdiately yield 

P ° Ut ~ { C' ext ) dP ™ + c C> xt c 
w\ S C' ext * C> 



c 

3 



Eq. (46) is covariant form of Eq. (11). 



(45) 



EKi W - -) (46) 
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4.4. Momentum-four vector for outgoing radiation - thermal emission 

If the particle is not irradiated, then, in accordance with Eq. (40), we have 
[(dp out )f = |f ^ - E - ^) I dr , (47) 

for outgoing radiation due to the thermal emission. 

5. Radiation: Energy-mass relation 

Now, we are interested in the relations between the corresponding part of radiation, 
incoming or outgoing, and it's mass. 

5.1. Incoming radiation 

Making an invariant of Eq. (42) 

(c dM ln f = dp? n d Pln „ , (48) 

one immediately obtains, on the basis of Eqs. (25) and (33) - b\ M = - 

dM in = ; (49) 

a summation over repeated upper and lower indices is always implied, e. g., dp^ n dpi n M 
= dp° in dp m o + dpj n d Pm i + dp 2 in dp ln 2 + dp\ n dp in 3 = {dpl n ) 2 - (dp} n ) 2 - (dp 2 m ) 2 
(dpi) 2 - 

On the basis of Eqs. (42) and (49), we can summarize: 

C Wi 

d Pin = ^^^dr, 

C Wl 

dM in = . (50) 

Eq. (49) states that energy of the incoming/incident radiation is nonzero, while (the 
invariant) mass of the radiation is zero. The zeroness of the mass of the incident radiation 
is understandable, since the radiation is produced by the parallel flux of photons moving 
in the same direction and orientation (see Okun 1989a, 1989b; Schroder 1990, p. 108: 
"For the invariant description of the inertial behaviour of a particle, only its rest mass 
m can be used, because it has the same value in all reference systems."). 

5.2. Outgoing radiation 

On the basis of Eqs. (25), (31)-(33) and (44), we can write 

3 C' 



dPout 
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u„ u; = c, j = 1, 2, 3 , 

„, ».M _ „2 



6i ju 



i,j = 1 to 3 



(51) 



where (5^ is Kronecker delta (Sij = 1 if i = j, 5^ = if i ^ j). 

Making an invariant of Eq. (44) or of the first equation in Eqs. (51), 

(c dM out f = dp out u dp^ ut , 

one immediately obtains, on the basis of Eqs. (51) and (52), 
_ w{ S C' ex 



(52) 



dM n 



ext 



w\ s c> x 



(53) 



On the basis of Eqs. (44) and (53), we can summarize: 



dE out — X c 



_j_ _ \ " / Cprj 
™1 h V 



X 



-1 



J = l 



dr , 



dM n 



X 



X 



r 

w\ S C' ea 



C e xt 



+ x-^ F> 



3 / n, 

E ( # + x- 1 Kj 



3 = 1 



dr 



(54) 



Eq. (54) states that energy of the outgoing radiation is nonzero, and also (the invariant) 
mass of the radiation is nonzero. The procedure of the calculation of the mass is consistent 
with the procedure discussed by Okun (1989a, 1989b). 
Eqs. (51) and (54) yield 

dE out = -f out (dM ou t) c 2 , 

7 OMt (dM out ) v out , 



dpout 

Vout 
^iout 

zc 



3 n, 

E(# + x- 1 n. 



3 = 1 




(zc)- 1 



zc 

~ZM ' 



= — v(S 



3 = 1 



'-"ext 



-7 



ZM = 



L VyL 



w\ S C' ea 



+ F> 



3 

E 



C 



-3-*- + X- 1 F' . 



(55) 



We have obtained the standard formulae: energy equals 7 x mass x c 2 , momentum 
equals 7 x mass x velocity. 
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5.3. Outgoing radiation - thermal emission 

Making an invariant of Eq. (47) 



[(c dM out ) e f = [(dp out ) e ]^ [(dpout) 



eJ ' 



(56) 



one easily obtains 



{dM out ) e 



\ 



El 



E 



dr . 



Eqs. (47) and (57) yield 



(dE out ) e = 
(dM out ) e = 



7 > dr , 



i 



3 /pi n 2 



Wj = 7 



where also Eq. (25) was used. 
Eqs. (47) and (58) yield 

(dE out ) e = le (dM out ) e c 2 , 
(dpout) e = le (dM out ) e v 
E' v 



e > 
3 



— 7 - 

c c 



7e 



(ZM) e 



(ZM), 



\ 



K) 2 - E( CF e, 



7:)} 



(57) 



(58) 



(59) 



We have obtained standard formulae: energy equals 7 x mass x c 2 , momentum equals 
7 x mass x velocity. 



6. Incoming and outgoing radiation - change of particle's mass 

As it was already pointed out just below Eq. (38), mass of the particle in the whole 
process of interaction with the radiation is conserved, i.e. dm/dr = 0; if the particle is 
not irradiated, then Eqs. (40)-(41) hold. Now, we are interested in the effect of incoming 
and outgoing radiation on the mass of the particle. 
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On the basis of Eq. (42), we can write for the change of four-momentum of the particle 
due to the incoming radiation 

W = Ml , 
dr dr 
dpi w\ S C' ea 



K , (60) 



dr c 
which leads to the result 

fdm_\ = wj S C' ext 1 
\dr] in c c ' 

Comparison of Eq. (61) with Eq. (50) yields 
dE^ = fdm\ c2 l_ 

dT \ dT J in W l 

and 

dp in ( dm\ ei 

c — . (63) 



6.2. Outgoing radiation 



On the basis of Eq. (51), we can write for the change of four-momentum of the particle 
due to the outgoing radiation 



dp* dp^ 



out 



dr dr 
dp» out w\SC' ext L *C> 



dr c 

3 



j=i \ ' 



- J2Ki -) • (64) 

This leads to the result 

fdm_\ = _ w\ S C' ext 1 _ 
V dr J t c c 

\ / out 

According to Eqs. (61) and (65): (dm/dT) out = — (dm/dr)^- Comparison of Eq. (65) 
with Eq. (54) yields 

dEoUt l 1 / ^Pr,j y-l JTl/ 

~*r ~ ~ i ^ - £ [c^ + x \ 

3 



dPout _ I ei | C'prj y _i rpi 



y ext 




X = r^U^Hii . ( 66 ) 



_ ™? 5 C' ext 
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6.3. Thermal emission 



As for the thermal emission alone, the relevant result is given by Eqs. (40)-(41): 

dm\ E' Q 
dr J e c 2 

Eqs. (40)-(41) immediately yield for the energy and momentum of the outgoing radiation: 



(67) 



Eq. (68) yields the standard formulae: energy equals 7 x mass x c 2 , momentum equals 
7 x mass x velocity in the case F' j = 0, j = 1 to 3. 



<iE / dm 



„2 



i=i 

3 



7. Discussion 

7.2. Physics of p^ n 

Eq. (11) is physically different from Eq. (11) in Klacka (2004) and also from Eq. (122) in 
Klacka (1992). In order to understand the physics behind the formulae, we will consider 
a simple case of spherically symmetric particle fulfilling C' pr 2 = C' pr 3 = C' abs = F' e l = 

K.2 = K, 3 = 0. now - 

For the sake of brevity, we will use dimensionless efficiency factors Q' x instead of cross 
sections C' x : C' x — Q' x A', where A' is geometrical cross section of the spherical particle, 
and, we will not use primed quantities. 

Let a beam of homogeneous incoming photons, interacting with the spherical particle, 
consists of N photons each of which carries energy e. The total incoming momentum is 
N (e/c) e and the elastically scattered photons are characterized by the total outgoing 
momentum (e/c) Y^j=\ e j' momentum of an individual photon is (e/c) ej, j = 1 to N. 
The momentum of the particle is given by conservation of the total momentum: N (e/c) 
Q pr s e = N (e/c) e - (e/c) ^=1 e h according to Klacka (1992, 2004). This yields 
Qprs e = e — X^jLi e j I N ■ Taking into account that the last term is an average which 
can be written as < cos 9 > e, we finally obtain Q pr ,i = 1 — < cos 9 >. However, this 
is not correct. We know that Q pr ,i A = C w s = C ext — < cos 9 > C sca and Q ext A = 
Qabs A + Q sca A = Q sca A, since we deal with the case Q a bs = 0, in our example. Thus, 
Q pr ,i = Qext ( 1 — < cos 9 > ) and the correct form of the conservation of the total 
momentum is N (e/c) Q pr ,\ e = Q ext [N (e/c) e — (e/c) X^Li e j ] ■■■■ (*)• This explains 
why the ratios of the pressure terms to the extinction term are present in Eq. (11), and, 
also, why the extinction cross section is present in energy in Eq. (11). The equation (*) 
is fully consistent with Eq. (11): p' t = ( C' ext /A' ) N' ( e'/c ) e', p' a = [1 - ( C' prA / C' ext 
)} p'i, where the dot denotes differentiation with respect to time. All is consistent with 
Eq. (5), too. 
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7.2. Center-of-momentum 

Short comment on Eqs. (55) and (59). We have defined velocities v out and v e for radiation 
in Eqs. (55) and (59). This in reality corresponds to the fact that we can introduce 
center-of-momentum frame, where the total momentum vanishes. Let v c be the velocity 
of the center-of-momcntum frame with respect to the reference frame in which an object 
(radiation, particle) of a mass dM moves with velocity v. Then 

alE = 7c (dM) c 2 , 
dp = 7c (dM) v c , 

= (69) 

and center- of -inertia of the system of the outgoing photons, defined as the point with 
radius vector Rs(t) — ( J2i E>i(t)ri(t) ) / ( J2i Ei(t) ), moves with the velocity v c given 
by Eq. (69) (see, e. g., Landau and Lifshitz 2005: §14, pp. 42-45, or, Fcrraro 2007: §6.5, 
pp. 146-148, and §6.7, pp. 153-155). 

7.3. Special case 

Now, we will deal with the special case C' pr2 = C pr3 = F' e - = for j = 1, 2, 3. This 
case corresponds to secular decrease of semi-major axis and eccentricity of neutral dust 
grain in the gravitational and electromagnetic fields of a central star, e. g., Sun, if the 
cross section C' prl is not a function of time (see Sees. 6.1 and 6.2 in Klacka 2004). Thus, 
we can call this effect as the (generalized) Poynting-Robertson effect. As in Sec. 6.1, we 
will use dimensionless efficiency factors Q' pr ^, Q' ext instead of cross sections C' pr l , C' ext : 
Q'pr.i = Cpr,i I ^ ■> Q'ext — C'ext I ^ i where A' is geometrical cross section of a sphere 
of volume equal to the volume of the particle. In reality, the case treated in this section 
holds only for particles with spherically symmetric distribution of mass. Dimensionless 
efficiency factors can be calculated on the basis of Mie theory (Mie 1908, Bohren and 
Huffman 1983). 

Equation of motion for the discussed case is given by reduction of Eq. (38): 
= (l/w 1 ,e 1 /w 1 ) , 

wi = 7 (1 — v ■ ei/c) , (70) 

where also Eqs. (18) and (33) were used. As a consequence, dm/ 'dr = (this corresponds 
to the condition dE' /dr = 0, where E' is energy of the particle measured in the rest 
frame of the particle). If the particle is not irradiated, then only thermal emission may 
exist and one has to use 



'd 

~J7 



, (71) 
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instead of Eq. (70) (sec Eq. 41). 

To first order in v/c, Eq. (70) yields 
d v S A' -. ( / v ■ e\ s 



dt mc ^ pr ' L l\ c 

It is worth mentioning to stress that the value Q' pr l does not depend on particle's orien- 
tation with respect to the incident radiation, but it's value may be time dependent, e. g., 
when optical properties change with the particle's distance from the central star (Klacka 
et al. 2007). 

The total process of interaction can be written in the form of the following equations 
(energies and momenta per unit time): 



E' = E' = S' A' 



'ext J 



P'o = (1 - 



' ext 



v\ = (E'Jc) e[ , (73) 

where the index "i" represents the incoming (incident) radiation, beam of photons, the 
index "o" represents the outgoing radiation. Eq. (73) is a fundamental condition for the 
validity of the Poynting-Robertson effect and it is a new condition. It differs from the 
condition p' a = ( 1 — Q' pr l ) consistent with the statements presented in Poynting 
(1903), Robertson (1937), Wyatt and Whipple (1950), Burns et al. (1979; Eqs. 10 and 
11) and in other papers. The condition represented by Eq. (73) does not exist in the 
literature and any statement about the outgoing radiation, published up to now, is not 
consistent with Eq. (73). 

Covariant form of Eq. (73) can be obtained from Eq. (46): 

dpLt _ ( 1 Q'prA Mi + W g A ' Q'ext Q'prS ^ ^ 

dr \ Q' ext ) dT c Q' ext c 

and incoming four-momentum of the radiation is given by Eq. (42) : 

d v - n = Wl S * ®' ext ftf dT . (75) 

Momentum-four vector for the outgoing radiation can be written in a short relativis- 
tically covariant form, on the basis of Eq. (44), as follows: 

d-Pout _ v I a 

~ " X Y 1 ~ Q'ext 

x ^ w\ S A' Q' ext ^ 
c 

vLt = (^f , Pout) • (76) 
and Eqs. (18), (32) and (33) can be used. 
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As for the incoming radiation, the relevant statements are represented by Eqs. (42) and 
(50), as for the four- momentum, and also, as for the energy and the mass of the radiation. 
Now, we will treat the outgoing radiation. Eq. (51) reduces to (see also Eq. 76) 

Q' V r.\ 



dPout 


wjS A' Q' ext 
c 




= c , 






hub? 


= . 



1 - 



Qext I Qe 



dr 



ext 



(77) 

Definition of the invariant mass represented by Eq. (52), dM out — dp out M dp^ ut / c 2 , 
leads to 

w\ S A' Q' ex 



dM 



Qpr.l 



Q'pr.l 



dr 



(78) 



° 2 ^ Q'ext \ Q'ext / 

if one uses results presented in Eq. (77); compare with Eq. (53). As a consequence of Eq. 
(78), one obtains a simple statement: 



<^i<2. 

Qext 

Eqs. (77) and (78) yield 
{ZC\ 



(79) 



dE out 
dpout 
(ZC\ 

(ZM\ 



(ZM\ 
1 

JZM\ 

u 



(dMout) c 2 , 

Qpr.l 



1 - 



Q'e 



ext 



^ + 

W l Q'ext 



pr.l V 

~~ 7 c 



(dM out ) c , 



? pr,l 



Q'e 



ext 



1 Q'pr,l 
™1 Q'ext 7 



Qpr.l 



\| Q'ex 

wi = 7 (l " 
see also Eq. (55). 



I 2 - 

v • a 



pr, . 
ext 



(80) 



v c = 



Q' 



As for the center-of-momentum frame, Eqs. (69) and (80) lead to 
Q' 7 v + (1 - Q') c ei/iui 



Q' 7 + (1 - Q') M 



(81) 



e:.r/ 



It can be verified that Eq. (81) yields 
v c = (l - Q') cei , if v = , 
t? = - (l - Q') c ei , if v c = ( 



v c = v , if Q' = 1 , 

Qpr,l 



Q' 



Q'e 



(82) 
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The case v = — ( 1 — Q' pr l / Q' ext ) c e i yields dE out = ( dM out ) c 2 , as can be verified 
on the basis of Eq. (80). 

If the particle is not irradiated, then the general Eqs. (58)-(59) reduce to 

(dE out ) e = 7 {dM out ) e c 2 , 
(dpout) e = 7 {dM out ) e v , 

(dM out ) e = Srdr . (83) 

Moreover, Eq. (59) or comparison of Eq. (83) with Eq. (69) immediately gives 

v c = v . (84) 

Eq. (84) states that the center-of-momentum frame for thermally emitted radiation is in 
the center of the spherical particle. 

7.3.2. Energy-mass relation for the particle 

At first, we will treat the incoming radiation. On the basis of Eqs. (61)-(62), we have 



This is not equivalent to the change of mass of the particle according to the "definition" 
c~ 2 dE in I dr ("mass equals energy divided by c 2 "), mainly if one takes into account 
that Eqs. (60)-(61) and (69) lead to v c = c e\: the relation between the momentum of 
the incoming radiation and the particle's mass is (see Eq. 63) 



We also remind of Eq. (49) for the mass of the incoming radiation: dM in / dr = 0. 
For the outgoing radiation, Eq. (66) reduces, for the spherically distributed mass, to 



This is not equivalent to the change of mass of the particle according to the "definition" 
j~ 1 c~ 2 dE out /dT. Changes of the particle's mass and outgoing momentum are related as 
follows: 




(85) 




(86) 




(87) 




(88) 
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Eqs. (87)- (88) yield, together with equation analogous to Eq. (69), that v c is given by 
Eq. (81), and, thus, we can write 



dr lc y \ «e j ^ ^ j 

Q' 7 v + (1 - Q') c ei/ioi 



«c = 



0' 7 + (1 - 0') M 

g/ ^ Qppl (89) 

The relations for E out and p OMt seem to be only a little modified standard relations 
(invariant multiplier \]Q'{2 — Q') is present). However, this is not true: the velocity is 
v c and it is not equal to the velocity of the particle v (except for the case Q' = 1): Eqs. 
(69) hold. Also, comparison of Eqs. (87)-(88) with Eq. (80) yields 

dM out 



dr 



y 



\ Q'eo 




(90) 



out 



Eqs. (69), (89)-(90) state that the formulation "mass equals energy divided by c 2 " holds 
for radiation (dE out , dM out ), but not for the mass of the particle (dE out , dm), if Q' ^ 1. 
It is worth mentioning that the special case Q' w ^ / Q' ext = 1, yielding p' Q = (generally 
believed to be the case of perfect absorption within geometrical optics approximation 
and which was treated by Robertson 1937, see also Robertson and Noonan 1968, p. 114), 
leads to very simple physics: v c = v, according to Eq. (89), dM out / dr — — (dm/dT) out , 
according to Eq. (90), and (dE) out — — 7 (dm) out c 2 . 

If the particle is not irradiated, then Eq. (41) holds: (dm/dr) e = — E' a / c 2 . 
Comparison with Eq. (83) yields: (dm/dr) e = — (dM out /dT) e . 



7.3.3. Mirror-like spherical surface and geometrical optics - specular reflection 

On the basis of Eq. (79) we know that < Q' w ^ j Q' ext < 2. Moreover, one could 
say that the special case of perfect absorption within geometrical optics approximation 
would yield Qp r ,i / Q'ext = 1; according to Eq. (73) and according to the conventional 
statement that the absorbed radiation is isotropically reemitted (in the particle's frame 
of reference) - p' a = 0. Moreover, a statement that Q' pr l / Q' ext — 2 should hold for 
spherical particle with mirror-like surface (totally reflecting sphere), can also be found. 
However, the statements presented in the last two sentences are not correct. In order to 
prove this conclusion, we will treat a nonrotating totally reflecting sphere, within a geo- 
metrical optics approximation. We will consider specular reflection of the electromagnetic 
radiation at the surface of the spherical particle, in the particle's rest frame of reference 
(for the sake of brevity, we will not use primed quantities, in this subsection). 

Let the center of the cartesian reference frame x — y — z be situated at the center of 
the spherical particle and the incoming beam of photons is characterized by unit vector 



22 J. Klacka: Electromagnetic radiation, motion of a particle 

of momentum p inc = — x = — ( 1, 0, 0). Any point on the surface of the particle, 
which interacts with the incoming radiation, can be characterized by coordinates of the 
surface position vector in the following way: 

R = Rn , 

n = (cos a cos 5, sin a cos<5, sin 6) , 

. 7T 7T. 

ae(--, + -), 

5 e (- \ + \) , (91) 

where R is particle's radius and n is unit vector. 

Now, we have to use the law of reflection: The angle of incidence equals the angle 
of reflection, and the incident and reflected rays are in the same plane. The angle of 
incidence <p is given by condition 

cos ip = I n ■ (— x) | = cos a cos 5 . (92) 

The plane of incidence is characterized by its normal unit vector N: 
(— x) x n 



N = 

\x x n\ 

(ft Qin A — Qin rv rriQ A I 

(93) 



\J sin 2 5 + sin 2 a cos 2 S 
The equation of the plane of incidence is given as N ■ R + k = 0, which yields k = — 
N ■ R = 0, according to Eqs. (91) and (93). The plane of incidence is identical to the 
plane of reflection: 

N ■ P out = , 

Pout ■ n = cos ip , (94) 

where also the condition " the angle of incidence equals the angle of reflection" was added 
(the case a = 0, S = yields p out = + x). Conditions formulated by Eq. (94) yield, on 
the basis of Eqs. (91)-(93): 

Pout = { x out} Vout} z out) J 

x out — (cos 2 a — sin 2 a) cos 2 S — sin 2 5 , 
cos a sin a cos 2 S 

Uout l.t %out) ~ 2 • 2~T 1 

sin a cos z o + sin o 
cos a cos S sin S 

Zout = (1 — Xout) —T~2 . 9 r • (95) 

sin a cos z o + sm d 
Unit vector along the small circle parallel with the equator e a , and, unit vector along 
the great circle normal to the equator e$, both localized at the position vector R, are 
given by formulae 



e a = (cos (a + , sin (a + , Oj 

= (— sin a, cos a, 0) , 
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eg = n x e a 

= (— cos a sin 8, — sin a sin 8, cos<5) , (96) 

since the relation for eg is a consequence of the relation e a x eg = n. Element of the 
length vector dl = dx x + dy y + dz z = R dn is, on the basis of Eq. (91): 

dl = R(— sin a cos 8 da — cos a sin 8 dS) x 
+ R (cos a cos 8 da — sin a sin 8 d8) y 

+ Rcos8 d8 z . (97) 
Components of the length vector along the vectors e Q and eg are: 
dl a = dl ■ e a — R cos 8 da , 

dig = dl e 5 = Rd8 . (98) 

Vector of the elementary/infinitesimal area is 

dA = (dl a e a) x (dig eg) — R 2 cos 8 da d8 n 

= R 2 dVLn , (99) 

where d£l is the solid angle element at the position (a, 8). The amount of energy (per 
unit time) hitting the area dA at the position R is S dA± , where S is the flux density of 
radiation energy (energy flow through unit area perpendicular to the ray of photons per 
unit time) and 

cL4^ = dA ■ x = (R 2 dfl) cos a cos 8 = 

= R 2 cos a cos 2 8 da dS , (100) 

if also Eq. (91) is used. The quantity dA± is projection of the area dA on the direction 
normal to x. 

On the basis of Eqs. (91) and (100) we can write for the incoming momentum per 
unit time 

dpi, 

— = - o j ^ Pl!> , 

r w/2 r ir/2 



- c S JdA ± Pi , 

— p inc / da d8 {R 2 cos a cos 2 8} 

c J-tt/2 J-tt/2 



rr/2 J-tt/2 

- c *r 2 p inc . (101) 



The outgoing momentum per unit time is, on the basis of Eqs. (95) and (100), 

1 V [dA ±Pout = 0. (102) 



dt c 

The mirror-like spherical particle obtains the following momentum per unit time: 



dp _ dp in _ dp ou t _ 5 2 
dt ~ dt dt ~ » n ^ mc ' 



tt R 2 p inc , (103) 
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which immediately follows from Eqs. (101)-(102). Comparison with Eq. (70) or (72) 
yields Q' pr l = 1. This is consistent with the results presented in van dc Hulst (van do 
Hulst 1981, p. 161: see Table 13 "Efficiency Factors for Totally Reflecting Spheres"). 
However, comparison of Eqs. (73) and (102) would yield Q' pr ,\IQ'ext = 1j an d, thus, 
Q' ext = 1. But this is not consistent with the results of the Mie theory presented in 
van dc Hulst (1981, p. 161 - Table 13). This seems surprisingly, since van de Hulst 
states on p. 223: "a smooth totally reflecting sphere with radius large compared to the 
wavelength scatters light by reflection isotropically" . Where's the problem? In order to 
understand physics of the results, we have to realize that the result presented by Eq. (102) 
is based on geometrical optics and the result does not respect all physical interactions, 
only reflection is considered. In reality, "extinction paradox" exists (van de Hulst 1981, p. 
107) and diffracted light plays a non-negligible role in treating the incoming and outgoing 
radiation, even for large particles (in comparison with a wavelength of the interacting 
light); "diffraction = small-angle scattering" (van de Hulst 1981, p. 107). We know that 
"the diffracted light gives a zero contribution" to the radiation pressure of large spheres 
(van de Hulst 1981, p. 225), but diffraction cannot be neglected in a separate treatement 
of the incoming and outgoing radiation. In order to be consistent with equations of the 
type of Eq. (52), i. e., Eqs. (53) and (78)-(79), we have to take into account diffraction, 
also. Thus, we have Q' w ^ = 1 and Q' ext = 2, both for perfectly absorbing sphere {Q' sca 
= 1 due to the diffracted light and Q' abs = 1) and for totally reflecting sphere (Q' sca = 1 
+ 1, 1 due to the diffracted light + 1 due to the reflected light, Q' abs = 0): 



Qpr,l _ 1 



(104) 



Q'ext 2 ' 

The result Q' pr l / Q' ex t = 1/2 holds for spherical particles not only with perfect ab- 
sorption, but, also with mirror-like surface characterized by specular reflection. Eq. (73) 
holds for the incoming and outgoing radiation. Eq. (78) yields 

dM out = V3 W * S 2 A ' dr . (105) 

One would await, on the basis of a simple idea about the isotropic light reemission or 
scattering (diffraction is neglected), the result without the numerical factor \/3 - the 
same result would be obtained on the basis of the following frequent idea (see also Sec. 
8): "Photons at rest do not exist! It is therefore somehow artificial to speak of the photon 
rest mass. It is more logical to use the relation m = h v / c 2 for the definition of the 
photon mass." (Demtroder 2006, p. 92). 

7.3.4. Summary 

We have already mentioned "the extinction paradox" (van de Hulst 1981, p. 107). It 
states that "a large particle removes from the incident beam exactly twice the amount 
of light it can intercept" (van de Hulst 1981, p. 107). As van de Hulst explains, the 
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explanation is in the fact that "(a) all scattered light, including that at small angles, is 
counted as removed from the beam, and (6) that the observation is made at a very great 
distance, i. e., far beyond the zone where a shadow can be distinguished". 

Let us take into account the values presented in Table 13 in van de Hulst (1981, p. 161). 
One can immediately realize that the values Q pr greater than 2 would be in contradiction 
with equations analogous to Eqs. (78)-(79), if p' = ( 1 — Q' pr ) p- would be the condition 
for the Poynting- Robertson effect. Thus, the case of "perfectly absorbing case" does not 
correspond to p' Q = within geometrical optics approximation plus diffraction (relevant 
physical processes), i. e., isotropic reradiation/reemission. Instead of this we have to 
use Eq. (11) and the "perfectly absorbing case" yields p' Q = ( 1 — Q' pr / Q' ext ) p- = 
(1 / 2) p- for large particles [2 tt x radius of the particle / (wavelength of the incident 
radiation) 3> 1]. The term "isotropic reemission" does not hold, strictly speaking, neither 
for small particle, nor for large sphere. Similar situation holds for totally reflecting sphere. 
[Although, within geometrical optics approximation, the term "isotropic reradiation" is 
used, which is based on the fact that "for large spheres we may separate the effects of 
geometrical optics (reflection and refraction) and the effect of diffraction" (van dc Hulst 
1981, p. 225).] 

We have found that the condition for the validity of the Poynting-Robertson effect 
is given by Eq. (73). The case of "perfectly absorbing spherical particle" is given by Eq. 
(104), or, by the condition p' a = 0.5 p- for large spherical particle. This condition differs 
from the condition used by Robertson (1937): it is generally accepted that the Robertson's 
case corresponds to p' a = 0. E. g.: "Consider a spherical particle with proper mass m and 
world- velocity u M which scatters (whether by reflection or by absorption and re-radiation) 
electromagnetic radiation isotropically in all directions relative to its proper system. ... 
The proper system x a of the particle at a given event is used for writing down the 
equation of motion. The force on the particle is then due only to the incident radiation, 
because the re-emitted radiation is isotropic." (Robertson and Noonan 1968, p. 114). 
Similarily, in his original paper, Robertson writes on p. 427: "the particle is assumed to 
possess spherical symmetry and to radiate isotropically in the rest system" of the particle. 
On the basis of this assumption, the author deals with geometrical optics approximation 
and he writes equation of motion in the form corresponding to our equation dp^ / dr = 
d Ptn / dr - dp% ut I dr, where dp^ n / dr is given by Eq. (75) and dp% ut / dr by Eq. (76). 
However, Robertson asserts that Q' prl — 1 and Q' ext = 1 (see Eqs. 2.7 and 2.8 on p. 427 
in Robertson 1937, and, also the statement that "the first term [(wfSA'/c) b%] represents 
the total power-force four-vector due to the incident radiation, with which it coincides 
completely in direction, and the second [{w\SA' /c) u^/c] that due to the radiation re- 
emitted at the rate [wfSA'c] "). The statements of Robertson are not correct, since Q' prl 
= 1 and Q' ext = 2, in reality (see Eq. 104). The Robertson's approach would yield, on 
the basis of equations similar to Eqs. (77)-(78) (with the transformation Q prl / Q' ext 
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— ► Q' prl , see e. g., Eq. 122 in Klacka 1992), that < Q' prl < 2, which would be an 
inconsistency between the Mie theory and the relativity theory. These considerations of 
Robertson are equivalent to that of Poynting (1903, p. 539: "... it is receiving a stream 
of momentum w\SA' /c directed from the sun. Its own radiation outwards being equal 
in all directions has zero resultant pressure." ... we have used symbols of our paper). In 
reality, the case p' Q = 0, considered by Poynting (1903), Robertson (1937), Wyatt and 
Whipple (1950, p. 558: "the process of absorption and re-emission produces no net force 
on a particle when one chooses to work with a stationary frame referred to the particle" 
... the part "absorption and" should be omitted in the formulation) and others, does not 
hold and the correct form is p' Q = 0.5 p\. Thus, we have understand the real physics of 
the P-R effect for the first time, now! We are aware that the effect of diffraction must 
not be neglected, even in the case of geometrical optics approximation. 

Finally, let us consider a large plane mirror of geometrical area A' characterized by 
specular reflection (totally reflecting plane) and normal incident radiation. We have C' abs 
= °> C 'sca = 2 A ', C 'ext = 2A ' and (cos0') = 0, C' prl = 2 A'. We can write 
P'o = ( i- c prA I C' ext ) p> = ... (**), 

and Eqs. (11) and (**) yield 

E[ = S> C' xt = 2 S> A', p'^iE'/c) e[, 

E' a = E[ = 2 S' A', p' a = 0. 

Eqs. (50) and (53) yield, then: 

dM m = 0, dM out = (2 S' A' I c 2 ) dr. 

On the basis of geometrical optics considerations without the effect of diffraction (C' pr l 
= 2 A', C' ext = A' would have to be used in Eqs. 50 and 54), one would await that p' Q 
= — Pi and ( dM out ) na ive = ( dM in ) na ive = - the idea that photon's mass equals 
its energy divided by c 2 would yield ( dM in ) ph = ( dM out ) ph = ( S' A' / c 2 ) dr. The 
result for dM out can be easily physically understood: the situation is similar to finding 
a mass of two photons, each of which travels in the same direction but orientations of 
their motions are different - reflected and diffracted photons are important, in our case. 

7.4. Orbital evolution 

In general, equation of motion for arbitrarily shaped particle is given by Eqs. (38) or 
(39), or, by Eqs. (40) or (41). Laboratory experiments carried out by Krauss and Wurm 
(2004) proved that nonspherical particles behave in a different way than spherical ones. 

The Poynting-Robertson effect yields the following property, as for the particle's or- 
bital evolution in gravitational and electromagnetic fields of a central star: semi- major 
axis and eccentricity secularly decrease (see Klacka 2004: Sees. 6.1 and 6.2 - mainly 6.2.3, 
for more details), if the optical properties of the particle do not change (an increase of 
the orbital elements can occur, if (3 increases: see Eq. 12 in Klacka 1993). This is fulfilled 
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under assumptions presented in Sec. 7.3, i.e. C' pr 2 = C' pr3 = F' e = and Eqs. (73) or 
(74) hold. The equation of motion is given by Eqs. (70) or (71), then. As a consequence, 
the equation of motion holds for particle with spherically symmetric mass distribution. 

As for presentations in current literature, the statements are not consistent with the 
discussed results. 



7.4.1. Accelerations of the spherical particle 

Relations for the accelerations of the spherical particle interacting with the incident 
electromagnetic radiation can be obtained from Eqs. (75)- (76). The relevant relations for 
the accelerations can be summarized in the following way. 

Eq. (75) yields for the four-acceleration of the spherical particle due to the incoming 
radiation: 



X s S ■"' «- , (106) 



m c 

or, to the first order mv/c 
fdv\ S A' Q' e 



\dtj in mc 
The change of the particle's mass is 
dm\ _w 2 1 SA> Q' ext 



dT , ,. (108) 

Eq. (76) enables to find the four-acceleration of the spherical particle due to the 
outgoing radiation. Eq. (76) yields for the force acting on the particle (dp^ /dr) ou t = — 
dPout/d-T, which yields for the four-accclcration of the particle 

du*\ ( _ Q^ rA \ / 



X = W ' - A ' ®' ext , (109) 
m c 

or, to the first order mv/c 
dv\ SA'Q' ext f i _ Q prA \ 



dt J out mC V Q ' ext J 

,{{1-^)^1}. 

The change of the particle's mass is 



dm\ w\ S A' Q' e 



ext 



dT / out 

Total four-acceleration of the spherical particle is 
fdu^\ fdu^\ (du»\ 



(111) 
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and the mass of the particle is conserved: {dm/ dT) part i C ie = {dm/dT) in + (dm/dT) out — 
0. To the first order in v/c 

If Q' prA /Q' ext = 1, then (dv/dt) out = (see Eq. 110). Further, < Q' prA /Q' ext < 2, 
according to Eq. (79). If Q' pr \/Q' ex t < 1, then the re-radiated light yields (dv/dt) out oc 
+ v. If Qpri/Q'ext > 1; then the re-radiated light yields (dv/dt) ou t oc — v. The case 
(dv / dt) ou t oc — v leads to a decrease of particle's total energy (kinetic plus potential) 
and to an increase of the orbital speed \v\, if the optical properties of the particle do not 
significantly change during the particle's motion around a star. 



7.4.2. Karttunen et al. (2007) 

Karttunen et al. (2007, p. 201) state: "When a small body absorbs and emits radiation, 
it loses its orbital angular momentum and the body spirals to the Sun." 

The statement holds only for spherical particle orbiting a star (Sun), under the as- 
sumption that the particle is characterized with spherically symmetric mass distribution 
and constant optical properties. Moreover, even when spherical particle does not absorb 
radiation, it may (secularly) lose orbital angular momentum (see also Sec. 7.3.3 as a 
very special case, or any elastic scattering). If optical parameters of the particle change, 
then the particle may not lose its (secular) orbital angular momentum. If the particle 
is nonspherical, there may not exist systematic secular decrease of particle's semi-major 
axis and eccentricity, and, moreover, particle's orbital plane may change. 



7.4.3. Carroll and Ostlie (2007) 

Carroll and Ostlie (2007, p. 806) describe the Poynting-Robertson effect in this way: 
"When particles absorb sunlight, they must re-radiate that energy again if they are to 
remain in thermal equilibrium. The original light was emitted from the Sun isotropically, 
but in the Sun's rest frame the re-radiated light is concentrated in the direction of motion 
of the particle. Since the re-radiated light carries away momentum as well as energy, the 
particle slows down and its orbit decay." 

The statement does not hold for arbitrarily shaped particles. It holds, partially, only 
for dust grains with spherically symmetric mass distribution and constant optical proper- 
ties. If such type of spherical particles is considered, then Eq. (74) holds for the outgoing 
radiation. Eq. (74) shows that "re-radiation" (but also scattering, in general) occurs both 
in directions dpi n /dt oc e\ and v. We can write, on the basis of Eqs. (77) and (78) (to 
the first order in v/c): dp out /dt = ( S A' Q' ext / c) x [ ( 1 — Q' ) ( 1 — v • ei / c) e.\ 
+ Q' v / c ] = ( S A' Q' ext I c ) x [ ( 1 - v ■ e x / c ) e x - Q' a + Q' ( ( v ■ d / c ) 
ei + v I c ) ], Q' = Q' pr l I Q' ext - Only the direction v is considered in the explanation 
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of Carroll and Ostlie (" ... in the Sun's rest frame the re-radiated light is concentrated 
in the direction of motion of the particle."). 

The statement that " the re-radiated light carries away momentum as well as energy" 
is correct. Mathematical description of this physical process is given by Eq. (76). But 
the complete sentence "Since the re-radiated light carries away momentum as well as 
energy, the particle slows down and its orbit decay." is physically incorrect, in general. 
Really, Eq. (76) yields for the force acting on the particle {dp* 1 /dT) out = — dpg Ut /dr, 
which yields for the four-acceleration of the particle results given by Eqs. (109)-(110). If 
Q'pr.i/Q'ext = 1) tnen (dv/dt) out = (see Eq. 110) and the particle does not slow down. If 
Q' pr ,i/Q'ext < 1' tncn tne re-radiated light yields (dv/dt) out oc + v. If Q' prt i/Q' ext > 1; then 
the re-radiated light yields (dv/dt) out oc — v. The case of "perfectly absorbing sphere", 
treated by Robertson (1937), corresponds to Q' pr ^/Q' ex t = 1/2 within geometrical optics 
approximation, and, thus, the re-radiated light yields (dv/dt) out oc + v, according to Eq. 
(110). 

7.4.4. Lissauer and Murray (2007) 

Lissauer and Murray (2007, p. 806) explain: 

"A small particle in orbit around the Sun absorbs solar radiation and reradiates the 
energy isotropically in its own frame. The particle thereby preferentially radiates (and 
loses momentum) in the forward direction in the inertial frame of the Sun. This leads 
to a decrease in the particle's energy and angular momentum and causes dust in bound 
orbits to spiral sunward. This effect is called the Poynting-Robertson drag. 

The net force on a rapidly rotating dust grains is given by 
Frad « [(£ Q P r A) / '(4 7t c r 2 )] [ (1 - 2 v r / c) r - (vq/c) © ] . (50-LM) 
The first term in Eq. (50-LM) is that due to radiation pressure and the second and the 
third terms (those involving the velocity of the particle) represent the Poynting-Robertson 
drag. 

From this discussion, it is clear that small-sized dust grains in the interplanetary 
medium are removed: (sub)-micronsized grains are blown out of the solar system, whereas 
larger particles spiral inward toward the Sun." 

(Explanation to symbols in Eq. 50-LM are given also on p. 806 of the paper: "A is 
the particle's geometric cross section, L is the solar luminosity, c is the speed of light, r is 
the heliocentric distance, and Q pr is the radiation pressure coefficient, which is equal to 
unity for a perfectly absorbing particle and is of order unity unless the particle is small 
compared to the wavelength of the radiation.") 

Now, we will treat the real physics. 

i) 

"A small particle ... ." (Lissauer and Murray 2007, p. 806) 
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The equation of motion corresponding to the Poynting-Robertson effect holds for 
(nonrotating) spherical particles. 

ii) 

"A small particle ... reradiates the energy isotropically in its own frame." (Lissauer and 
Murray 2007, p. 806) 

General assumption/condition for the validity of the Poynting-Robertson effect is for- 
mulated by Eq. (73). This equation shows that p' Q = if and only if Q' i/Q' ext = 1. 
Thus, the statement of the authors, presented above, corresponds only to a special com- 
bination of optical properties of the spherical particle and the wavelength of the incident 
electromagnetic radiation. The authors probably assume perfectly absorbing spherical 
particles within geometrical optics approximation (moreover, diffraction is ignored by 
the authors), but this case is characterized by the condition Q' pr i/Q' ext = 1/2, as it is 
given by Eq. (104), and p' Q = (1/2) p\. 

iii) 

"The particle therefore preferentially radiates in the forward direction in the inertial 
frame of the Sun." (Lissauer and Murray 2007, p. 806) 

According to Eq. (76), we can write for the outgoing radiation for the P-R effect (first 
order in v/c is considered): 

d Pout /dt = (SA> Q' ext /c ) {(1 - Q' pr JQ' ext )(l - v ■ ei /c) ei + (Q' prA /Q' ext ) v/c] 

Since Q' pr _\/Q'ext 7^ 1 in general, also radial term is present in dp out /dt, not only 
" forward direction" , in the inertial frame of the Sun. 

As for the perfectly absorbing spherical particle: the standard approach Q' pr l = Q' ext 
= 1 would yield dp out /dt = ( S A'/c ) v/c = ( 5" A'/c ) {{v R /c)e R + (v T /c)e T }, while 
the correct result is Q' prl = 1, Q' ext = 2, and, dp out /dt = ( S A'/c ) {e R + (v T /c)e T }; 
the velocity vector v is decomposed into radial component v R e R and transversal compo- 
nent VT&T-, ejj = ei. 

iv) 

" The particle thereby preferentially radiates in the forward direction in the inertial frame 
of the Sun. This causes dust in bound orbits to spiral sunward." (Lissauer and Murray 
2007, p. 806) 

As it is evident from Eq. (110), the special case Q' pry i/Q' ext — 1 yields (dv/dt) out — 
and the reradiated energy does not cause dust to spiral sunward. Dealing with more 
realistic case Q pr ^/Q' ext 1) one obtains: 

a) if Q' pr _il Q'ext < 1) t nen the re-radiated light yields (dv/dt) out oc + v, 

b) if Q' pr _ilQ'ext > 1) tnen t ne re-radiated light yields (dv/dt) out oc — v. 

See also Sec. 7.4. f. Moreover, if optical properties of the dust change, then the P-R effect 
can cause spiralling outward from the Sun. 
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v) 

"The net force on a rapidly rotating dust grains is given by ... (Eq. 50-LM)" (Lissauer 
and Murray 2007, p. 806) 

In reality, Eq. (50-LM) holds for nonrotating spherical particle, as it is evident from 
derivations presented in our paper (see also Sec. 7.3.3). If we would like to consider 
rapidly rotating arbitrarily shaped dust particle, we have to consider general equation 
of motion (see Eq. 39) and the corresponding averaging over rotational motion has to 
be done (see, e. g., Klacka and Kocifaj 2001). Krauss and Wurm (2004) have given an 
experimental evidence that arbitrarily shaped and rapidly rotating dust grain moves in 
a different way than it corresponds to the Poynting- Robertson effect. 

vi) 

Equation (50-LM) contains general factor Q pr (= Q' pr 1 ). However, this is in contradiction 
with the statement " A small particle in orbit around the Sun absorbs solar radiation and 
reradiates the energy isotropically in its own frame." (Lissauer and Murray 2007, p. 806). 
The isotropically reradiated energy within geometrical optics approximation (diffraction 
is not considered in this formulation) corresponds to the condition Q' pr i/Q' ext — 1/2 
(see Eq. 104 in Sec. 7.3.3), in reality. Thus, if the authors state that the essence of 
the Poynting- Robertson effect is "particle reradiates the energy isotropically in its own 
frame", then they have to use Q pr = 1. 

It is also important to stress that the process "particle reradiates the energy isotrop- 
ically in its own frame" (Lissauer and Murray 2007, p. 806) is only an approximation 
within geometrical optics and that, in reality, the fundamental condition for the case is 
given by the equation p' a = (1/2) p\ and not p' a = 0. The last equation considers the effect 
of diffraction (calculation of mass of the outgoing radiation within relativity approach 
requires consideration of the effect of diffraction), while the phrase "particle reradiates 
the energy isotropically in its own frame" does not consider the effect of diffraction. 

vii) 

" The first term in Eq. (50-LM) is that due to radiation pressure and the second and the 
third terms (those involving the velocity of the particle) represent the Poynting- Robertson 
drag." (Lissauer and Murray 2007, p. 806) 

The effect of electromagnetic radiation on spherical dust particle is given by Eq. (70). 
It corresponds to radiation pressure. Eq. (70) contains all the terms present in Eq. (50- 
LM). Thus, Eq. (50-LM) corresponds to the radiation pressure - not only its first term, 
as the authors state. Physics does not allow to divide Eq. (70) into several parts which 
can be treated as separate physical phenomena. 

viii) "From this discussion, it is clear that small-sized dust grains in the interplanetary 
medium are removed: (sub)-micronsized grains are blown out of the solar system, whereas 
larger particles spiral inward toward the Sun." (Lissauer and Murray 2007, p. 806) 
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The conclusion made by the authors does not follow from the discussion presented by 
the authors. Even if we take into account all the corrections presented above and that the 
correct and fundamental result is given by Eq. (50-LM). The authors' conclusion may be 
violated if one takes into account that optical properties of the particles may change. 

7.4.5. Aberration (of light) 

It is often stated that the aberration of the incoming light is responsible for the term 
— v/c in the Poynting-Robertson effect (see, e. g., Dohnanyi 1978, p. 562; Leinert and 
Gran 1990, p. 226; Gran 2007, p. 632). Is it true? If the incoming radiation would be 
responsible for the term — v/c in the P-R effect ("radiation falls preferentially on the 
leading edge of the orbiting particle and acts as a drag force" - Festou et al. 2004, p. 729), 
then this term should be present in pi = dp in /dT. However, Eqs. (42) and (75) show that 
dpi n /d,T is proportional to 6i = e\/w\. Thus, the term — v/c does not come from the 
incoming radiation. In reality, the term — v/c is present in dp out /dT (Eq. 76). But it's 
presence is not due to the aberration of light. It's existence is caused by conservation of 
particle's mass, when considering both the incoming and the outgoing radiation. Really, 
if we would require, e. g., dE'/dr = {w\ S Q' pr A' / c) a' and dp' /dr = {w\ S Q' pr A' 
/ c) e[, then Eq. (14) would yield dp/dr = {w\ S Q' pr A' / c) [ &i + ( a' - 1 ) u/c ], 
where bi — e\/w\, u = 7 v, wi — 7 ( 1 — v ■ e\ / c ), and, a 1 is a constant. The term 
( a' — 1 ) "fv/c does not correspond to the aberration of light. As a similar example, we 
can mention Eq. (40). 

7.5. Analogy between mechanics and electromagnetism? 

We have discussed the effect of electromagnetic radiation (photons) on perfectly absorb- 
ing and totally reflecting spherical and planar particles. We have stressed that even within 
the geometrical optics approximation the effect of diffraction is important, in the process 
of interaction between the particle and the electromagnetic radiation. As a consequence, 
the relation between the outgoing and incoming momenta per unit time (Eq. 73) reduces 
to p' — 0.5 p\ for the above mentioned spherical particles (Eq. 104). The numerical 
factor 0.5 considers also the effect of diffraction which is conventionally neglected (e. g., 
Lissauer and Murray 2007). 

Thus, we know that "for large spheres we may separate the effects of geometrical 
optics (reflection and refraction) and the effect of diffraction" (van de Hulst 1981, p. 
225). But even in the case of geometrical optics approach, the simultaneous action of 
these effects cannot be neglected, if one wants to be consistent with the relativity theory. 
If we take into account a plane particle (arbitrarily shaped particle, in general), all these 
effects are also important. 
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The effect of reflection of the electromagnetic radiation on a planar object is often 
explained in an analogy to mechanical process. As an example, we can mention two 
well-known textbooks on fundamental physics. Halliday et al. (2008, p. 900) write about 
the comparison between the radiation effect on perfectly absorbing and totally reflecting 
object: "Instead of being absorbed, the radiation can be reflected by the object: that is, 
the radiation can be sent off in a new direction as if it bounced off the object. If the radi- 
ation is entirely reflected back along its original path, the magnitude of the momentum 
change of the object is twice that given above ... (the case for total absorption)" . "In the 
same way, an object undergoes twice as much momentum change when a perfectly elastic 
tennis ball is bounced from it as when it is struck by a perfectly inelastic ball (a lump of 
wet putty, say) of the same mass and velocity." Similarly, Jewett and Serway (2008, p. 
963) write: "Electromagnetic waves transport linear momentum as well as energy. ... In 
this discussion, let's assume the electromagnetic wave strikes the surface at normal inci- 
dence ..." "Momentum transported to the perfectly absorbing surface has a magnitude 
p." The authors continue (Jewett and Serway 2008, p. 964): "If the surface is a perfect 
reflector (such as a mirror) and incidence is normal, the momentum transported to the 
surface ... is twice that given by (the magnitude p). That is, the momentum transferred 
to the surface by the incoming light is p and that transferred to the surface by the re- 
flected light also is p v . In other words, the momentum transferred to the surface by the 
incoming light equals the momentum transferred by the reflected light. If the mechanical 
process would be a good analogue to the electromagnetic process, then the relation p' Q 
= ( 1 — C' prl I C' ext ) p- would reduce to p' a = — p-, for the perfect (both mechanical 
and electromagnetic) reflection. But the reality is different. The mechanical process is 
characterized by the values C' prl = 2 A' , C' ext = A' , while the electromagnetic process, 
including the effect of diffraction, yields cross sections C' prl = 2 A', C' ext — 2 A', where 
A' is geometrical cross section of the object. 

Summarization of the incident and outgoing momenta per unit time can be done in 
the following way (A' is geometrical cross section of the planar surface, u' and c are the 
speeds of the incident mechanical particle and photon with respect to the surface, and, 
5" is the rate of flow of energy, i. e., the rate at which the energy flows through a unit 
surface area perpendicular to the direction of the particle or photon propagation), for 
normal incident energy: 

i) perfectly absorbing planar surface: 

• mechanics: p \ = ( S' A'/c ) ( v! / c ) e[, p' a = 

• electromagnetism: p[ = 2 ( 5" A'/c ) e[, p' = ( S' A'/c ) e[ 

[ 1 x ( 5" A'/c ) e[ comes from the effect of diffraction, simultaneously in p[ and p' a ] 

ii) perfectly reflecting planar surface: 

• mechanics: p\ = ( S' A'/c ) ( u' / c ) e[, p' a = - ( S' A'/c ) ( v! / c ) e\ 

• electromagnetism: p\ = 2 ( S' A'/c ) e[, p' = 
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[ as for the p' Q : the effect of diffraction ( S' A'/c ) e[ is compensated by the effect of 
reflection - ( 5" A'/c ) e[ ] . 

Finally, we will present a short explanation to the results of mechanics. Four-momentum 
of the incident (classical) mechanical particle is p% P = ( E' MP /c, p' MP ). Four- momentum 
per unit time is p'f = ( E'Jc, p\ ) = n' u' A' p% p = ( n' u' A' E' MP /c ) x ( 1, u'/c ), 
since p' MP = E' MP u' / c 2 , u' = u' e[. Defining 5" = n' u' E' MP , we obtain p'/ 1 = ( 5" 
A' I c ) x ( 1, v! e[ j c ). 

8. On the equivalence principle of mass and energy 

One of the well-known results of the special relativity theory is the equivalence of mass 
and energy (more properly: the relation between mass and energy) . This " universal equiv- 
alence principle of mass and energy" is often formulated in the following way (or a little 
different, but equivalent form): "there does not exist a mass without energy (and vice 
versa) and every change in energy is connected with a corresponding change of inertial 
mass" (Schroder 1990, p. 114). It is stated that it holds also for a photon (light) and 
an inertial mass of the photon equals m p hoton — E/c 2 , where E is energy of the pho- 
ton. A simple "proof can be found in Kittel et al. (1962; Sec. 12.6, Eqs. 58-59, p. 400): 
P = "niphotonC and E/c = p, which yields E — m v h oton c 2 , p is momentum of the photon. 
(Serway et al. 2005, p. 95: "The photon has zero mass, but its effective inertial mass m,i 
may reasonably be taken to be the mass equivalent of the photon energy E, or m, = 
E/c 2 = hf jc 2 . The same result is obtained if we divide the photon momentum by the 
photon speed c: nij = p/c = hf jc 2 ? Analogously, Gasiorowicz 1974, p. 466 states: "A 
photon of energy E has gravitational mass E/c 2 .") The principle is usually proved by 
the thought experiment, coming back to Einstein (see, e. g., Einstein 1906, p. 633; Beiser 
1969, Sec. 2.6; Bernstein et al. 2000, pp. 81-82; Chow 2008, pp. 272-273; compare also 
Okun 1989b, p. 34). The idea is following: "Classical electromagnetism tells us that light 
carries both momentum and energy. Because of this momentum, light exerts a pressure. 
We can think of the radiation as being composed of massless particles - any object that 

moves with the speed of light must have no mass - carrying momentum A box on 

a frictionless surface has a radiation emitter on one end and an absorbing surface on 
the other. When a burst of radiation is emitted, the box recoils, stopping only when the 
radiation is absorbed at the other end. The center of mass will not move if there is mass- 
energy equivalence." (Bernstein et al. 2000, p. 81). As a result, the relation for a beam 
of parallelly spreading photons is presented: E — m c 2 , where E is energy of the photons 
(electromagnetic radiation), c is the speed of light and m is the mass corresponding to 
the energy E (the result of the thought experiment is: " The energy of the radiation is 
equivalent to a mass according to the Einstein formula: E — m c 2 ." Bernstein et al. 2000, 
p. 82). 
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However, the following question arises: What is the physical sense of the mass for 
a beam of masless photons? According to Schroder (1990, p. 106): "The inertia of the 
particle is described by an invariant quantity, the mass m of the particle, more properly 
called its rest mass". On the other hand, the author writes: "According to the universal 
equivalence principle of mass and energy, there does not exist a mass without energy 
(or vice versa) and every change in energy is connected with a corresponding change of 
inertial mass." (Schroder 1990, p. 114). Thus, an inconsistency emerges: the inertia of 
a beam of massless particles, photons, is described by an invariant quantity, the mass 
of the beam m(, eam = 0, and, the beam possesses and energy Eb eam which should be 
equivalent to a mass according to the Einstein formula mb eam — Eb eam / c 2 . 

8.1. Thought experiment 

In trying to better understand the situation, we wil present detailed calculations for a 
thought experiment. The thought experiment will be similar to that discussed above, but 
instead of a box of mass M we will consider two surfaces of masses M/2 on a frictionlcss 
surface. The surfaces are penpendicular to the frictionlcss surface. They are at rest in 
an inertial frame of reference K, at the beginning, with initial positions xl = —L/2 
and x R = +L/2. The right surface emits a particle with energy E (measured in the 
reference frame K) towards the left surface, i. e. direction and orientation of the motion 
of the particle is characterized by the unit vector — x. Let the emission of the energy E 
corresponds to a mass decrease m R of the right surface. The absorption of the energy E 
produces a mass increase mt of the left surface. 



8.1.1. Emission of a massless particle - particle with zero mass 

After emission of a massless particle (photon) by the right surface and its absorption by 

the left surface, the center of inertia of the system does not change (see, e. g., Eq. 14.6 
in Landau and Lifshitz 2005): 

E L r L + E R r R = 0. (114) 
We can rewrite Eq. (114) in detail: 

1L + m L ^j + v L t^j = 7fl - m R ^j + v R ^ + v R tj , (115) 

where t = holds for the moment of the photon absorption by the left surface, 7l = 
1 / \f\ — (vl/c) 2 , j r = 1/ yjl — (v R /c) 2 ; the multiplicative factors c 2 present in energies 
E L = 1l{M/2 + m L ) c 2 and E R = j R (M/2 - m R ) c 2 are omitted. Eq. (109) yields 

7,(f +rn L )= lR (^ -m R ) (l + 2 V -f) (116) 

for t = 0. Since the center of inertia does not change also for t > 0, Eq. (115) yields 

(M \ (M \ , N 

7i t + m L\ v L = 7fl v - rn R \ v R . (117) 
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Eq. (117) corresponds to the laws of conservation of momentum for the moments of 

emission and absorption: 

( M \ E 

Tii h m R) V R = — > ( 118 ) 



( M \ E 

1L [ y + m L \ v L = — . (119) 

The laws of conservation of energy for the emission and the absorption are of the 

following forms in the inertial frame of reference K: 

(M \ , M , . . 

1R f y - rn R \ c 2 + E=— c 2 , (120) 

(M \ , M , , N 

7l (^y + mi J c = y c + £ • ( 121 ) 

Eqs. (118) and (120) yield 




or, 



(l - =«) ^ . ,123, 



During the process of the emission of the photon(s), the mass rriR (more correctly: the 
energy ttir c 2 ) has been changed into the energy E and kinetic energy of the surface of 
mass M/2 - m R [kinetic energy = (7^ - 1) (M/2 - m R ) c 2 ]; m R w [ 1 + ( S/c 2 ) / M] 
-E/c 2 . This is in agreement with Eq. (120) and this is the physical interpretation of the 
Eqs. (118), (120) and (122)-(123). Mass of the photon(s) with the energy E is zero. 
Eqs. (119) and (121) yield 




™l = ^t Wl + 2 -L— - 1 , (124) 



or, 



E 



= (l + ^) m L c 2 . (125) 



During the process of the absorption of the photon(s), the energy E has been changed 
into the mass (more correctly: rest energy ttll c 2 ) and motion of the surface of mass 
M/2 + m L [kinetic energy = (j L - 1) (M/2 + m L ) c 2 ]; m L w[l-( E/c 2 ) / M } E/c 2 . 
This is in agreement with Eq. (121) and this is the physical interpretation of Eqs. (119), 
(121) and (124)-(125). 



8.1.2. Summary 

We can make the following statements: 

i) Photon(s) (light) with the energy E transferred mass mi. 
i) Mass of the photon(s) (light) is zero. 
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Proof of the statements: 

i) Initial mass of the left surface was M/2, it's mass after absorbing the photon(s) 
is yj '{energy / c 2 ) 2 — {momentum / c) 2 = ^/[7l(M/2 + mi,)] 2 — [jl{M/2 + m^v^fc] 2 = 
M/2 + rriL- The mass tul is given by Eq. (117), or, approximately, tol w [ 1 — {E/c 2 ) 
I M] E/c 2 . 

ii) Mass of the photon(s) equals {energy /c 2 ) 2 — {momentum / c) 2 = y 7 {E/c 2 ) 2 — {p/c) 2 



= ^/{E/c 2 ) 2 ~[{E/c){~x)/c] 2 = ^J{E/c 2 ) 2 - {E/c 2 ) 2 = 0, where also the fact that the 
photon(s) moves with the speed c in the direction and orientation — x, was used. 

The statements are consistent with the statements presented in Taylor and Wheeler 
(1992, p. 228): "Energy without mass: photon", "photon moves with zero mass". 

8.1.3. Emission of a particle with non-zero mass 

After emission of a particle with non-zero mass by the right surface and its absorption 
by the left surface, the center of inertia of the system does not change, see Eq. (114). We 
can rewrite Eq. (114) in detail: 

1L + m^j + v L tj = 7fl - m R ^j + v R ^ + v R tj , (126) 

where t = holds for the moment of the absorption of the particle by the left surface, 
7l = 1/ \A — ( v l/c) 2 , "fR = 1/ \/l — {vr/c) 2 , u is the speed of the particle in the frame 
of reference K; the multiplicative factors c 2 present in energies El — Jl{M/2 + tul) c 2 
and E R = j R (M/2 - m R ) c 2 are omitted. Eq. (126) yields 

„(f - m „) ( 1 + 2 2f) 

for t = 0. Since the center of inertia does not change also for t > 0, Eq. (126) yields 

(M \ (M \ 
7i t + m L) v L ="fR\ — m R \ v R . (128) 



2 J ' \2 

Eq. (128) corresponds to the laws of conservation of momentum for the moments of 

emission and absorption: 

(M \ r s 

lR \ ~2 _ mR I VR = 7 " m U ' *- 129 ^ 

1l ( y + m L j v L = 7« m u , (130) 



where m is mass of the particle and 7 U = 1/y 1 — {u/c) 2 . 

The laws of conservation of energy for the emission and the absorption are of the 
following forms in the inertial frame of reference K: 

( M \ i i M 2 . , 

7« y - ™fl c + lu m c = — c , (131) 

( M \ 9 M , , , 

1L [ y + m L j c 2 = — c 2 + lu mc 2 . (132) 
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Eqs. (129) and (131) yield 




2 E ' C2 



M/2 \M/2 



(133) 



or, 



n 777 r> 777 o , . 

m « c ^T7 c ' m < m « ' ( 134 ) 



(l - — ) m.c 2 + -mc 

"m 

During the process of the emission of the massive particle, the energy wir c 2 has been 
changed into the energy E and kinetic energy of the surface of mass M/2 — m,R [kinetic 
energy = (7^ — 1) (M/2 — mn) c ]. This is in agreement with Eq. (131) and this is the 
physical interpretation of Eqs. (129), (131) and (133)-(134). Mass of the massive particle 
with the energy E is m. 

Eqs. (130) and (132) yield 

E = 7„ m c 2 , (135) 
or, 



E = V + n) mL c - M mc 



2 , m i 771,2 2 „ 
= mi, c + — c , m < mL . (1j6) 

During the process of the absorption of the massive particle, the energy E has been 
changed into the mass ttll (more correctly: rest energy (?) and motion of the surface 
of mass M/2 + rriL [kinetic energy = (7^ — 1) (M/2 + tol) c 2 ]. This is in agreement with 
Eq. (130) and this is the physical interpretation of the Eqs. (130), (132) and (135)-(136). 



8.2. Summary 

The important statements, as for massive particles, can be summarized as follows: "The 
theory of relativity leads to the important conclusion that the energy of a particle (with 
m 7^ 0) at rest is mc 2 . The relations Eq — mc 2 , E = mc 2 / 'y/l — (u/c) 2 for a particle 
at rest and a particle in motion, respectively, are the famous Einstein formulae. They 
express the equivalence of mass and energy." (Schroder 1990, p. 113). As for the massless 
particle, it moves with the speed c and, thus, it cannot be at rest with respect to any 
frame of reference: rest energy of the particle is E n = 0. The inertia of the particle is 
described by an invariant quantity, the mass of the particle. The mass of a photon is 
zero. 
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The term "inertial mass of a photon" is of no physical sense: the speed of the photon is 
constant and the photon cannot be accelerated or decelerated, in a vacuum. Similarily, the 
term "gravitational mass of the photon" is of no physical sense - no (thought) experiment 
exists which would prove the correctness of the gravitational character of the photon. (See 
also Ugarov 1969, p. 248). 

Finally, we will reproduce several sentences on the meaning of E = m c 2 from Sachs 
(2007, p. 84): "It has been said by many physicists and philosophers that the formula 
E = m c 2 [we have added the subscript 0] means that 'mass is equivalent to energy'. 
This is philosophically false. It is not what Einstein said when he derived this relation. 
What he said was that 'the inertial mass of matter is a measure of its energy content'. In 
physics, as Newton originally postulated, the inertial mass of matter is, by definition, a 
measure of its resistance to a change of its state of rest or constant motion. The energy 
of matter, on the other hand, is by definition, the capability of this matter to do work. 
Thus, mass and energy are totally different concepts! What should be said, instead of 
saying that mass is equivalent to energy, is that mass (the inertia of matter) is a measure 
of the capability of this matter to do work (its intrinsic energy)." 

Results obtained in Sees. 4-7 are consistent with the previous statements (see also 
Eq. 105 with comment). 

8.3. Comparison of the results 

We want to compare the results obtained in Sec. 8.1 with the results of Sees. 2-7. 
8.3.1. Relativity and optics from Sees. 2-7 

Let us consider an incoming electromagnetic radiation which interacts with arbitrarily 
shaped particle. The incoming radiation is characterized by its energy and momentum. 
The relevant relations are given by Eqs. (42), (50), (60)-(63), or, by Eqs. (85)-(86) for 
spherical particles. In order to be the results easily comparable to the considerations 
presented in Sec. 8.1.1, we will discuss a totally absorbing planar surface of geometrical 
cross section A' and normal incident radiation, within geometrical optics approximation. 
On the basis of Sees. 7.3.3 (Eqs. 104-105 and the text above and below them) and 7.5, 
we can immediately write 



C' 



ext 



2 A' 



c: 



sea 



A' , 



dMi. 




= 2 S' A' , 



(137) 
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where S' is the flux of radiation energy and four-momentum per unit time is pf = ( 
E'Jc, p- ); C' ext and C' sca are cross sections for extinction and scattering. 

If the absorption and thermal emission are present, then the outgoing radiation is 
characterized by the following relations 



K 


= El 




Cpr, 1 


= A' , 




p'o 


■0 


C'ext J 


dM out 


= V3 


S' A' , 

n 2 dT ' 



(138) 



if the results from Eq. (137) are used. 

Now, let the thermal emission does not exist (thought experiment). Then 

E'o = 0, 
P'o = , 

C'pr, 1 = ^ A' , 

dM out = , (139) 

since p'f p' a ^>0,p'» = ( E'Jc, p' ),p' = (l- C p7 . 1 / C' ext ) p\. Eqs. (50) and (61)-(63) 
yield 



dM in = 0, 
dm \ w\ S C' ext 
dr 



tn 



dEi n ( dm\ 2 1 
c 



dr \dTj m wi 
dp in ( dm\ ei 



wi = 7„ (1 - u/c) , 
C' ext = 2 A' , (140) 

where also some results of Eq. (137) have been added. 
8.3.2. Relativity from Sec. 8.1.1 

Sec. 8.1.1 presents result for an increase of mass of a planar absorbing surface if it's mass 
is M/2 and the increase of it's mass is generated by an absorption of a photon of energy 
E with respect to the surface. In order to be able to compare the results of Sec. 8.1.1 
with the results presented in Sec. 8.3.1, we have to generalize the considerations of Sec. 
8.1.1 for the case of a moving absorbing surface. 

Let us consider a photon of energy E in the laboratory frame (inertial frame of refer- 
ence K) in which the absorbing surface is moving with the speed v, in the same direction 
and orientation as the photon is travelling (also negative orientation may be considered: 
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it is sufficient to put — v into the formulae presented below, in this subsection). We are 

interested in the increase of mass of the absorbing surface. The law of conservation of 

momentum for the absorption is of the following form in the inertial frame of reference 

K (compare with Eq. 119): 

/ M \ E M , N 

7i V + m L ) v L = h J v — v . (141) 



Similarily, the law of conservation of energy yields (compare with Eq. 121): 

1L f y + m L J c 2 = E + 7 „ — c 2 . (142) 

Eqs. (141)-(142) yield 



M \ E/c 2 



2 1 y M/2 j 

wi = 7„ (l - . (143) 

During the process of the absorption of the photon(s), the energy E has been changed 
into the mass tul (more correctly: rest energy tul c 2 ) and motion of it, and, part of the 
energy E was changed into the change of the kinetic energy of the surface of mass M/2. 
This can be written as: E = m L c 2 + (il - 1) (M/2 + m L ) c 2 - (j v - 1) (M/2) c 2 , or, 
E = 7l m L c 2 + (7 L - 7„) (M/2) c 2 . 

The limiting case M — > oo of Eq. (143) yields 
E 



lim rriL = Wi 

M— >oo 



(i - i) . 



Eq. (144) can be written as 
E' 

lim m L = — , (145) 

M^oo c 2 

where E' is energy of the photon in the frame of reference of the absorbing surface. 



8.3.3. Relativity from Sec. 8.1.3 

Sec. 8.1.3 presents result for an increase of mass of a planar absorbing surface if it's mass 
is M/2 and the increase of it's mass is generated by an absorption of a massive particle 
of energy E = -f u m c 2 with respect to the surface. In order to be able to compare the 
results of Sec. 8.1.3 with with the results presented in Sec. 8.3.2, we have to generalize 
the considerations of Sec. 8.1.3 for the case of a moving absorbing surface. 

Let us consider a massive particle of mass m and energy E in the laboratory frame 
(inertial frame of reference K) in which the absorbing surface is moving with the speed 
v, in the same direction and orientation as the massive particle is travelling, u > v (also 
negative orientation may be considered: it is sufficient to put — v into the formulae 
presented below, in this subsection: since u > also u > — v). We are interested in the 
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increase of mass of the absorbing surface. The law of conservation of momentum for the 
absorption gets the following form in the inertial frame of reference K (compare with 
Eq. 130): 

f M \ M . . 

1l I y + m L I v L =-/ u mu + 7„ — v . (146) 

Similarily, the law of conservation of energy yields (compare with Eq. 132): 

(M \ , , M , . . 

7l ( y + m L 1 c 2 = 7„ m c 2 l]„yc 2 . (147) 

Eqs. (146)-(147) yield 

= t |v + 2 ^ I 1 - -) m + \m) ~ l ) ' (148) 



Eq. (148) may be rewritten also in the following forms: 



E/c 2 


( m 


M/2 + 


\M/2 


m 


/ m 


M/2 + 


\M/2 



2 



m L = ^- { dl + 2w u + [ — — ) - 1 



w u = 7 „ (l - ^) , 
£ = 7„ m c 2 , 
7«' = 7« 7m ^1 - — J 



«' = — • ( 149 ) 

1 - U V J c l 

The limiting case M —> oo of Eqs. (148)-(149) yields 
lim wil = 7„/ m 

M— >oo 

E 

= W u -J , 

<r 

- _ " ^ \ 

7m' — 7m 7m ^t c 2 y ' 

Wm = 7m (l - > 

E = 7„ m c 2 . (150) 

Eq. (150) can be written as 
E 1 

lim m L — — , (151) 

M^oo C 2 

where i?' is energy of the massive particle in the frame of reference of the absorbing 
surface. 



8.3.4. Comparison of Sees. 8.3.1 and 8.3.2 

Let us consider a continuous beam of photons, each of which is bearing the energy E 
in the laboratory frame K. Let the absorbing surface is moving with the speed v in the 
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laboratory frame, in the same direction and orientation as the photons are travelling. 
The increase of mass of the absorbing surface, per unit time, is 

f dm \ S' C' ext , 1co , 

UO. =mL x -tf*- (152) 

\ / in 

Now, we have to use Eqs. (19), (23) and (144). Eq. (152) leads to 



dm\ E w 2 S C' 

■ = wi — x 



dr J .„ c 2 wi E 

_ S C' ext 



(153) 



The result obtained in Eq. (153) is identical to the result presented in Eq. (140) (see also 
Eq. 61). Thus, we have found a consistency between considerations presented in Sees. 2-7 
and Sec. 8.1.1. The term w\ present in Eqs. (143)-(144) is due to the Doppler effect; the 
same holds for the relations dE in /dT and dp in /dT in Eq. (140). 

8.3.5. Comparison of Sees. 8.3.2 and 8.3.3 

Absorption of a photon bearing the energy E in the laboratory frame K leads to the 
increase of mass of the absorbing surface moving with the speed v in the laboratory frame. 
The results are given by Eqs. (143)-(145). Similarily, absorption of a massive particle (the 
particle with a non-zero mass) bearing the energy E in the laboratory frame K also leads 
to the increase of mass of the absorbing surface moving with the speed v in the laboratory 
frame. The results are given by Eqs. (149)-(151). Eqs. (145) and (151) are equivalent. 
Moreover, comparison of Eq. (144) with Eq. (145) yields 

E' = 7, (l - I) E , (154) 

for the photon, and, comparison of Eq. (150) with Eq. (151) leads to 

E' = l v (l - ^) E, (155) 

for the massive particle moving with the speed u in the laboratory frame of reference K. 
In both cases, E is the energy of the incident photon / massive particle in the laboratory 
frame of reference and E' is the energy of the photon / massive particle measured in 
the reference frame of the absorbing surface. Eq. (154), which holds for the photon, 
corresponds to the Doppler effect. Analogous equation for the massive particle is Eq. 
(155) and it, formally, reduces to Eq. (154) in the limiting case u — > c. 

According to Eq. (150), we have limM^oo«i,L = 7«' m, or, limM^ooWL = w u E / 
c 2 . The quantity m is the mass of the massive particle and it is a relativistic invariant 
(VPm I c )- The quantity E / c 2 is not equal to m: m < liniM^oo < E / c 2 . Thus, 
it is of no sense to use E / c 2 as a mass of the particle. As a consequence, also E / c 2 in 
Eq. (144) does not represent a mass of a photon. These statements are consistent with 
Sec. 8.2 
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9. Cross sections for arbitrarily shaped particles 

We have already discussed, in Sees. 5, 6, 7.3 and 8, the relations between four-momentum 
of the radiation and the mass of the radiation or the change of the mass of the objects 
interacting with the radiation. Application of these physical considerations to astro- 
physics was shown mainly in Sec. 7.4 for the case of dust grains with spherically sym- 
metric mass distribution. Discussion in Sec. 7.4 explained the physical access to the 
Poynting- Robertson effect, which is standardly considered in evolution of interplanetary 
dust grains. This section will discuss the relation between four-momentum of the outgoing 
radiation and its mass for arbitrarily shaped dust grains. 

Let us consider that an arbitrarily shaped dust particle is irradiated by the incoming 
radiation of a source of radiation. The relevant relations for the four-momentum of the 
outgoing radiation and its mass are given by equations presented in Sec. 5.2. In what 
follows, we want to obtain more information on cross sections (for pressure components, 
extinction, scattering) than can be obtained from observations and measurements. 

As follows from Eqs. (53)- (54), the following inequality holds: 

X = Wl S °' ext . (156) 
c 

If we neglect the thermal emission force F' e , then Eq. (156) reduces to 
f" 3 / f" \ 2 

^ext j—i \ ^ext ) 

Current situation is that we do not know simultaneous values of the cross sections C' j 
(j = 1, 2, 3), C' ext , C' sca and C' abs from observations or measurements (Krauss and Wurm 
2004). Moreover, if the size of the particle is larger than about 5 micrometers, then it is 
not possible to make numerical calculations for the above presented cross sections of the 
particle. Thus, there do not exist methods which determine, simultaneously, the values of 
cross sections C' pr ^ (j = 1, 2, 3), C' ext , C' sca and C' abs , at present. However, we will show 
that some information on the cross sections can be obtained on the basis of equations 
presented in Sees. 2, 3.4 and 5.2, or, Eq. (157). 

Let us suppose that measurements can yield ratios of the cross sections for radiation 
pressure, i. e., the ratios C' j / C' prl (j — 2, 3) are in disposal. It is useful to rewrite 
Eq. (157) into the form 

o < ^r 1 < 3—^ 2 • ( 158 ) 

Eq. (158) reduces to Eq. (79) for the case C' pr 2 = C pr 3 = 0. On the basis of Eqs. (9)-(10) 
and (38), we can write C' pr l = C' ext — {< cos 9' > C' sca ) (the outer symbol for the mean 
denotes weighting over stellar spectrum and it corresponds to the "bar" symbol used 
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above the letters, standardly used in this paper). Since < cos 9' > < 1, we immediately 
obtain 

^F- > Sf 1 - 1 • ( 159 ) 

The equation C' ext = C' abs + C' sca can be rewritten into the form C' ext / C' prA = C' abs / 
C' prl + C' sca I C' pr l , which, together with Eq. (159) yields 

C' 

< < 1 . (160) 

pr,l 

Similarily, the relation C' ext / C' prA = C' abs / C' pr ^ + C' sca / C prA and Eq. (160) yield 
< Sf^ < ^££1 . (161) 
Eqs. (158)-(161) offer some information on the cross sections of extinction and scattering. 



10. Summary and conclusions 

The paper derives and presents relativistically covariant equation of motion for dust 
particle under the action of electromagnetic radiation - see Eqs. (38) and (40). It yields, 
as special cases, the results obtained by Einstein (1905) and Robertson (1937). As for 
most frequent applications to systems in the universe (e. g., meteoroids in the Solar 
System, dust particles in circumstellar disks), equation of motion in the form of Eq. (39) 
is sufficient. 

The general equation of motion reduces to the Poynting-Robertson effect for irra- 
diated particle with spherically symmetric mass distribution. This yields linear relation 
between the incoming and outgoing momenta (per unit time) of the radiation, in the 
proper reference frame of the particle: p' Q — ( 1 — Q' pr l / Q' ext ) Pi, where Q' pr l and 
Q' ext are dimensionless efficiency factors for the radial direction radiation pressure and 
extinction, integrated over the stellar spectrum (see Eq. 73, or Eq. 74 for covariant for- 
mulation). A simple condition < Q' pr l / Q' ex t < 2 is obtained from the invariant mass of 
the outgoing radiation (see Eqs. 77-79). The case of "perfectly absorbing" spherical dust 
particle corresponds to Q' pr ^i / Q' ex t = 1/2, if radius of the particle is much larger than 
the wavelcngth(s) of the incident radiation: the effect of diffraction cannot be neglected 
within geometrical optics approximation. While the condition p' Q = ( 1 — Q pr l / Q' ext 
) p'j is consistent with the condition dM out = y 'dp out M dp^ ut / c 2 , these fundamental 
relations are not consistent with the explanations of the Poynting-Robertson effect in the 
literature. 

The Poynting-Robertson effect is generated by simultaneous action of the Doppler 
effect, change of concentration of photons (the corresponding terms w\ x w\ are present 
at flux density of radiation energy S in equation of motion of spherical particle, in Eq. 
70) and the aberration of light, together with the laws of relativity theory. However, the 
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term — v/c docs not correspond to the aberration of light. It's existence is caused by 
conservation of mass of the particle. The term 1/wi in the four- vector 6^ (see Eq. 70) 
comes from the Doppler effect (see Eqs. 17-18) or from the change of concentration of 
photons (see Eqs. 21-22): Lorentz transformation of the zcro-th component of a four- 
vector is relevant, in relation between e! Y and e\. 

Covariant four-accelerations for the incoming and outgoing radiation are given by Eqs. 
(106) and (109). Moreover, the outgoing radiation accelerates the particle if Q' pr 1 /Q' ext 
< 1 (see Eq. 110), i. e., also the perfectly absorbing sphere treated by Poynting (1903), 
Robertson (1937) and others. 

As for the incoming radiation, the formulae for the energy and the mass of the ra- 
diation are given in Eq. (50). The statement that the mass corresponds to the energy 
through the relation "mass equals energy / c 2 " is not correct, because the photon is 
massless (see also Sec. 8). The relation between the incoming energy and the increase of 
the particle's mass is presented in Eq. (62) (Eq. 63 presents the result for momentum). As 
a consequence, the formula "energy equals 7 x mass x c 2 " does not hold. The Einstein's 
principle of equivalence of incrtial mass and rest-energy holds (Einstein 1999, p. 43). 

As for the outgoing radiation: if the particle is under the action of incident radiation, 
then the formulae for the energy and mass of the radiation are given by Eq. (54), see 
also Eq. (55): the simple formula "energy equals 7 times mass times c 2 " holds (see also 
Eq. 55). The relation between the outgoing energy and decrease of the particle's mass 
is presented in Eq. (66) for general case, and in Eq. (89) for particle with spherically 
symmetric mass distribution. As a consequence, the formula "energy equals 7 x mass x 
c 2 " does not hold, in general. The Einstein's principle of equivalence of inertial mass and 
rest-energy holds. 

Eqs. (40)-(41) present effect of thermal emission alone, if the particle is not irradiated. 
As a consequence, the isotropic thermal emission does not influence acceleration of the 
particle. Eq. (41) (Eq. 67) also shows that the formulation "energy equals mass times 
c 2 " holds between the decrease of particle's mass and the thermally emitted energy in 
the rest frame of the particle. General equation is given by Eq. (68) and the thermal 
force disturbs the saying "energy equals 7 x mass x c 2 ". The relation between energy 
and mass of the radiation is given by Eqs. (58) and (59): the standard simple formula 
"energy equals 7 x mass x c 2 " holds. 

On the basis of the covariant formulations presented in the paper, we were able to 
show that three of the current statements on the essence of the Poynting-Robertson effect 
are not correct. We have explained the important points in the statements. 

The relativistically covariant formulations presented in the paper enable to under- 
stand: (i) physics of the Poynting-Robertson effect (see also Sees. 7.3.4 and 7.4), which 
is frequently used in studies of orbital evolution of cosmic dust particles, and, (ii) why 
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the analogy on momentum transfer to the surface of an object in mechanics and electro- 
magnetism does not hold (see Sec. 7.5). 

We have also shown an application of four-momentum (and mass) of the outgoing 
radiation for arbitrarily shaped dust particle on obtaining some information on cross sec- 
tions, unknown from experiments, observations or theoretical/numerical solutions. The 
condition is presented in Eq. (158). As a result, the non-radial components of radiation 
pressure cross sections decrease the ratio C' x / C' ext in comparison with the spherical 
particles. A more general condition is given by Eq. (156). The condition can help in 
better understanding of optical properties of cosmic dust particles (special cases of the 
condition, applied to the spherical and planar particles, enabled us to understand physics 
of the interaction between the incident electromagnetic radiation and the particles, see 
Eqs. 79, 104, Sees. 7.3.4, 7.4, 7.5). 

The physics presented in the paper has direct implications for understanding of dis- 
tribution and evolution of cosmic dust grains in various astrophysical systems. 
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